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We study the evolution of the dynamical spin structure factor in a spin-1 antiferromagnetic Heisen- 
berg chain which is randomly doped with spin-1/2 moments. Using stochastic series expansion 
Quantum Monte Carlo simulations combined with the Maximum Entropy method, we monitor the 
crossover from the spin-1 chain, dominated by a single resonance at the Haldane energy, to the 
spin-1/2 chain with a continuum of spinon states which diverges at zero-frequency. Upon increasing 
the doping level, spectral weight is rapidly transferred from the Haldane peak to lower energies. 
If the exchange couplings between the spin-1/2 substituents and the spin-1 sites of the host are 
sufficiently small, finite-frequency bound states are observed below the spin gap for small doping 
concentrations. 
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It is well established that antiferromagnetic (AF) 
spin-1 Heisenberg chains have a low-energy spin gap 
A ~ 0.41J, leading to properties significantly different 
from quasi- long-range-ordered spin-1/2 AF Heisenberg 
chains.Eju The effects of non-magnetic impurities on such 
low-dimensional quantum spin liquids have been the fo- 
cus of considerable recent attention. In particular, it was 
shown that open ends in segmented Haldane chains in- 
duce low-energy states below the spin gapjj and that va- 
cancies in AF Heisenberg ladders liberate free spin de- 
grees of freedom in their immediate vicinityO'Q. At suffi- 
ciently low energies, RKKY-like interactions among these 
"pruned" spins lead to the formation of networks de- 
scribed by effective models of r|aBjdom spins interacting 
via random exchange couplings Jj'lI The effects of substi- 
tutions with non-vanishing magnetic moments tend to 
be more complex, depending on the spins of the dopants 
as well fSir^^i their exchange couplings with the host 
material.EijEJ We address this issue by studying the evo- 
lution of the dynamical spin structure factor of spin-1 
Heisenberg chains when doped with static spin-1/2 mo- 
ments, and by monitoring the crossover from a com- 
pletely gapped Haldane spectrum to that of a spin-1/2 
chain upon complete substitution. 

There have been several experimental and theoreti- 
cal studies of the effects of spin-1/2 impurities in the 
spin-1 chain. For example, the substitution of Cu^"*" for 
Ni2+ in the compound Ni(C2H8N2)N02(C104) (NENP) 
introduces spin-1/2 impurities on the original Ni sites, 
which have a woak ferromagnetic coupling to the S = 1 
chain-end spinsila On the other hand, substituting Ca^+ 
for Y'^+ in Y2BaNi05 introduces carriers in the chains 
which become localized on the oxygen sites between the 
Ni^+ ions. This breaks the superexchange path and re- 
places it biL,a direct exchange coupling to the spin-1/2 
impurities. Ej In this latter compound, inelastic neutron 
scattering experiments found a substantial increase of the 
spectral function below the Haldane gap, indicating-.the 
creation of states below the energy of the spin gap.Ej 

Numerical diagonalization and density matrix renor- 
malization group (DMRG) techniques were used to 



study the emergence of such impucitjiiinduced low-energy 
states in the Heisenberg modeliljO The main conclu- 
sion of these studies is that localized excited states ex- 
ist below the Haldane gap as long as the exchange cou- 
plings J' between the spin-1/2 and spin-1 sites are small 
compared to the couplings J between the original spin-1 
sites J, namely if J'LJ < 0.45. Other works have con- 
sidered mobile holesJiil'Ej In addition to enhanced low- 
energy spectral weight in the doped compound, these 
studies also observe low-temperature spin-glass behavior, 
due to impurity-induced frustrating ferromagnetic inter- 
chain interactions.Ej 



FIG. 1. Spin-1 AF Heisenberg chain with spin-1/2 mo- 
ments. Exchange couplings between spin-1 sites are indicated 
by solid lines, whereas the bonds between spin-1/2 and spin-1 
sites are denoted by dashed lines. 

He|Cfi__we concentrate on the static case discussed in 
Refs.t30 and apply stochastic, series expansion Monte 
Carlo (SSE QMC) simulatimistll in combination with the 
maximum entropy methodH to obtain the spectral prop- 
erties of Haldane chains randomly doped with spin-1/2 
sites. This approach allows us to study the evolution 
of the dynamical spin structure factor upon increasing 
the spin-1/2 concentration far beyond the limit of low 
impurity concentration. The doped spin-1 AF Heisen- 
berg chain pictured in Fig. 1 is modeled by a mixed-spin 
Heisenberg model. 
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H — JiSi • S 
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where Si denotes a spin-1 or spin-1/2 moment. In gen- 
eral, three inequivalent types of bonds exist in a mixed- 
spin chain: (i) bonds between two spin-1 sites, (ii) bonds 
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between two spin- 1/2 sites, and (iii) bonds between spin- 
1 and spin-1/2 sites. In the mixed spin model of Eq. 1, 
we use J as the original coupling of the Haldane chain 
with Ji — J if Si and Si+i are both spin-1, and oth- 
erwise Ji — J', such that the couplings between spin-1 
and spin-1/2 sites are J'. Since we are interested in the 
generic properties of this problem, in our simulations we 
concentrate on representative parameter sets. |— . 

The significance of J' has been emphasized in Ref.t2l 
for the case of isolated spin-1/2 impurities. Consider first 
the two extreme limits J'/J = and J'/J = oo. If the 
spin-1/2 dopants are only weakly coupled to the spin-1 
moments of the host, an effective Hamiltonian for the 
three-site cluster formed by the spin-1 chain-end spins 
and the impurity spin can be derived, correctly describing 
the emergence oLbound states below the Haldane gap of 
the spin-1 chain.E£l In the strong coupling limit, the cou- 
pling to the remainder of the chain destabilizes the local 
excitations. A critical couplpa-j'atio {J'/J)crit = 0.45 
separates these two regimes. llJO Therefore, in the fol- 
lowing we discuss separately the cases of weak coupling, 
which we illustrate with the parameter choice J' = O.IJ, 
and strong coupling, for which we take J' = J. 
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FIG. 2. Spectral function ^(Tr, uj) of the spin-1 chain doped 
with static spin-1/2 moments for various doping concentration 
p. The exchange couplings between spin-1 and spin-1/2 sites 
are J' = O.IJ (weak-coupling limit). 

In order to observe the evolution of the dynamical 
spectrum with doping cepcentration, we have performed 
SSE QMC simulationalJ and measured the imaginary- 
time spin-spin correlation function 
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considering the dominant AF wave vector kp^ tt by using 
the Monte Carlo estimator derived in Ref.EJ. We stud- 
ied chains of up to 500 sites for up to 500 realizations of 
the random distribution of spin-1/2 substituents down to 



temperatures T ^ O.OIJ. The maximum entropy analyt- 
ical continuation procedure was used to obtain the spec- 
tral function ^(Tr, oj) from the impurity-averaged data for 
the imaginary-time correlation functions.cS 

Let us first consider the case J'/J = 0.1, where the 
emergence of localized states below the spin gap is ex- 
pected fipm-DMRG calculations in the single impurity 
problemJlj'El In Fig. ^ we show the evolution of the 
spectral function S{tt,iv) upon increasing the concentra- 
tion p of spin-1/2 dopants. For the pure spin-1 chain, 
p = 0, a dominant resonance is observed at the Hal- 
dane gap energy w = A = 0.41J,|-eonsistent with ear- 
lier results using similar methods £3 Upon introducing 
spin-1/2 sites, the spectral weight of the Haldane peak is 
strongly suppressed and its width broadens. In addition, 
a significant amount of spectral weight appears at low 
but finite energies, indicating the emergence of impurity- 
induced states below the gap edge. In fact, the shape 
of the spectral density for p — 0.01 is consistent with 
Ref.t3, i.e. bound states with energies of approximately 
aJ'/2 and a3J'/2 appear in the low-energy spectrum, 
where a = 1.064 is a renormalization factor for the effec- 
tive jSpin- 1/2 chain-end spins Sl,r next to the impurity 
site.Ey These localized states are the lowest excited states 
of the effective Hamiltonian Heff =|-aJ'(SL • S' + S' • Su), 
where S' denotes the impurity spin.llj The maximum en- 
tropy method does not allow us to resolve the detailed 
structure of the spectrum, but clearly indicates the in- 
crease of spectral weight at this energy scale with in- 
creasing impurity concentration. Upon a further increase 
of the spin-1/2 concentration beyond p > 0.05, a zero- 
frequency peak emerges out of these low-energy contribu- 
tions. Simultaneously, thcJJaldane peak appears to move 
towards larger frequencie£j. The spectrum at larger im- 
purity concentrations {p > 0.2) exhibits features similar 
to those of the pure spin-1/2 chain, i.e., a divergence at 

= and no remaining peak structure at the energy of 
the Haldane gap. 
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FIG. 3. Same as Fig. 2, but with isotropic exchange cou- 
pling constants, J' = J, between all spin sites. 
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Turning to the isotropic case, J' = J, DMRG cal- 
culations for single spin- 1/2 impurities do not exhibit 
any low-energy states below the Haldane gap.ll3 Instead, 
an extended state with an energy of the order of A 
is reported. We performed SSE QMC simulations for 
the isotropic case and obtained the dynamical struc- 
ture factors shown in Fig. |^. Concentrating first on 
the low-energy properties of this spectral function, we 
can confirm that in the limit of small impurity concen- 
trations p < 0.06 no indications of low-lying impurity 
states are observed, consistent with the DMRG results. 
Upon further increasing the impurity concentration, a 
peak appears at zero energy, which develops into the 
zero-frequency singularity characteristic of the gapless 
spinon continuum of states in the pure spin- 1/2 Heisen- 
berg chain, i.e. S'(7r,w) oc Vhiuj/uj. Analogous to the 
discussion of the weak-coupling limit, spectral weight is 
shifted from the Haldane peak to lower energies. How- 
ever, in the isotropic case no bound states are observed 
at finite frequencies below the spin gap. Furthermore, 
the depletion of the Haldane resonance is less rapid than 
in the weak-coupling limit, indicating that the presence 
of weaker bonds has a more dramatic effect on the redis- 
tribution of spectral weight than the presence of spin-1/2 
substituents. 




FIG. 4. Spectral weight I of the Haldane peak, normalized 
to the spectral weight in the clean limit, lo- The circles repre- 
sent the isotropic case J' — J, whereas the squares denote the 
weak-coupling limit at J' = O.IJ. The dashed and dot-dashed 
lines are linear extrapolations. 



coupling. For example, at p = 0.5, the Haldane res- 
onance is still visible in Fig. || (J' = J), whereas in 
Fig. ^ (J' — O.IJ) this feature cannot be resolved be- 
yond p = 0.2. This critical value for J' = O.IJ is ob- 
tained by a linear extrapolation of / (dashed line shown 
in Fig. Q). In the case of isotropic couplings the spec- 
tral weight of the Haldane peak initially shows a rapid 
decrease, similar to the weak-coupling limit. However, 
beyond a doping concentration p ^ 0.06 at which zero- 
frequency singularities emerge in the dynamical structure 
factor (Fig. ||), / decays at a significantly slower rate. It 
is obviously more difficult to deplete the extended states 
at the energy of the Haldane gap for isotropic couplings 
than dissolving_the sub-gap bound states in the weak- 
coupling limit .E3 A linear extrapolation of / (dot-dashed 
line in Fig. ^) suggests a critical doping concentration 
p w 0.8 for isotropic couplings. Beyond this concentra- 
tion the dynamical structure factor resembles that of the 
pure spin-1/2 AF Heisenberg chain.tJ 

In conclusion, we have used a combination of stochastic 
series expansion Quantum Monte Carlo simulations and 
Maximum Entropy analytical continuation techniques to 
study the effects of random spin-1/2 substitutions on 
the dynamical response function of spin-1 AF Heisen- 
berg chains. In the limit of small doping concentrations, 
the dynamical spin structure factor is found to agree very 
well with results of earlier DMRG and exact diagonaliza- 
tion studies. If the AF couplings between the spin-1/2 
dopants and the spin-1 sites of the host are weak, bound 
states are formed below the Haldane gap A ~ 0.41J. 
On the other hand, in the case of isotropic coupling con- 
stants there are extended states at the energy of the spin 
gap. Upon increasing the doping concentration, a prolif- 
eration of low-energy states is observed at the expense of 
the spectral weight of the dominant resonance at w = A 
of the pure spin-1 AF Heisenberg chain. In the limit of 
weak coupling, the spectral weight of the Haldane peak 
disappears rapidly. However, for isotropic couplings the 
transfer of spectral weight to lower energies is more grad- 
ual. Here, the gapless spinon continuum spectrum of the 
pure spin-1/2 AF Heisenberg chain is recovered beyond 
a critical doping concentration p « 0.8. 
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In order to quantify this statement, the dependence 
of the spectral weight / at the Haldane energy is plot- 
ted in Fig. ^ as a function of the doping concentration. 
Here, a local Gaussian fit of the spectral function cen- 
tered around the peak position was used to estimate 
the spectral weight from the height and the width of 
the approximant. It is found that in the isotropic case 
remnants of the Haldane peak can be resolved at much 
higher spin-1/2 concentrations than in the limit of weak 
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